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Abstract. After a brief review of the Hartree-Fock-Bogoliubov method with Skyrme effec- 
tive interactions, we show how it can be applied to the description of various nuclear systems, 
from finite nuclei to neutron-star crusts. 



1 Introduction 

The global description of properties of finite nuclei over the entire nuclear chart 
requires theoretical methods as microscopic as possible and at the same time com- 
putationally tractable. It has long been recognised that self-consistent mean-field 
methods with effective nucleon-nucleon interactions can be very successfully ap- 
plied for this purpose (TJ. 

In particular, the Brussels-Montreal group has developed a series of nuclear 
mass models based on the Hartree-Fock-Bogoliubov (HFB) method with Skyrme 
effective forces [2]. The model parameters are fitted to essentially all the available 
atomic mass data, with the constraint to reproduce several nuclear-matter properties 
as obtained from microscopic calculations using realistic nucleon-nucleon forces. 
In our latest model HFB- 17 [3|, we have achieved our best fit ever to essentially 
all the available experimental data, the rms deviation for the set of 2149 measured 
masses of nuclei with N and Z > 8 [4| being only 0.581 MeV. Our model was 
also constrained to fit the the equation of state of neutron matter, as calculated by 
Friedman and Pandharipande [5 | for realistic two- and three-body forces. Besides 
our model reproduces the 1 So pairing gaps in both symmetric nuclear matter and 
neutron matter from the recent Brueckner calculations of Cao et al. |6|. Because of 
these additional constraints, our model can be used to reliably extrapolate beyond 
the neutron drip line and study astrophysical environments like for instance the inner 
crust of neutron stars . 



2 Skyrme-Hartree-Fock-Bogoliubov mass models 

All our HFB mass models are based on a conventional Skyrme force of the form 
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where r.y — Ti — rj, r = (r, + fj)/2, pij = — i7t(V, — Vj)/2 is the relative 
momentum, and Pg- is the two-body spin-exchange operator. Following the usual 
practice, we consider a different force in the pairing channel. The latter acts only 
between nucleons of the same charge state q (q — n or p for neutron or proton, 
respectively) and is given by 

vf T (n, rj ) = v*«[p n (r),p p (r)] S(r l3 ) , (2) 

where v* q [p n , p p ] is a functional of the nucleon densities. 

Assuming time-reversal invariance, the ground-state energy can be written as 
the integral of a purely local energy-density functional £hfb (r) which depends on 

(i) the nucleon density (denoting the spin states by a — ±1), 

p q ( r ) = Y A?( r ' cr ; r ' cr ) . ( 3 ) 

CT=±1 

(ii) the kinetic-energy density (in units of h 2 /2M q ), 

r q (r)= ]T f dV«5(r-r')V-VV,(r, ff ;r' l( T) 
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(iii) the spin-current density, 

Jq(r) = -i Y J d 3 r' S(r-r')Vp q (r,a-y,a') x ov ff 
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and (iv) the abnormal density, 



p«( r ) = X! / 5 <z( r ' cr ; r ' cr ) » ( 6 ) 

<r=±l 



where er^/ denotes the Pauli spin matrices. In turn the normal and abnormal density 
matrices, p(r, a; r' , a') and p(r, cr; r' , cr') respectively, can be expressed as ll9l[T0ll 
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where (r, a) and ^/4i ( r J <T ) me two components of the quasiparticle wave- 
function. Minimizing the HFB energy with respect to tp[f (r, <r) and 4>2i(r, a) un- 
der the constraints of fixed particle numbers leads to the HFB equations ll9lfT0l 

h' q {r) aa , A q {r)5 aa ,\ (^{r,a'j\ _ 
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where A g are Lagrange multipliers. The single particle Hamiltonian h' q (r) atT > and 
pairing field Z\ g (r) are given by 

h! q (r) a , a = -V • B q (r)V8 aa , + U q (r)S aa , - iW q (r) • V x tr a . a (10) 

and 

A q (r) = ~v*i[Pn(r),P P (r)Mr). (11) 
The single particle fields appearing in Eq. ( TTOb are defined by 

Bq( r) = d£ ™f^ , Uq (r) = V • ^ HFB(r) 



dr q {r) ' ^ ' 9p 9 (r) 0(Vp,(r)) ' 

Expressions for these fields can be found for instance in Ref. [8|. In the absence of 
pairing, the HFB equations (0 reduce to the Hartree-Fock equations. 

In homogeneous matter the HFB equations (0 can be readily solved. In partic- 
ular, the pairing field is uniform and obey the well-known BCS gap equations (see 
for instance Appendix B of Ref. (8l) 

A q = - ^Bl'\* « [p n , Pp ]A q [ de . ^ =, (13) 
8n Ja ^{e-\Y + Al 

where the subscript A is to indicate that the integral has to be regularized by intro- 
ducing a cutoff. In our latest models HFB-16 |8| and HFB-17 Q, we have inverted 
the gap equation ( fT3l ) at each neutron and proton density in order to deduce the 
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Figure 1. Differences between experimental and calculated masses as a function of the neu- 
tron number N for the HFB-17 mass model. 

effective pairing strength v 7Tq [p n , p p ] from a given gap function A q (p n , p p ). In par- 
ticular, for the model HFB-17 we have taken the 1 Sq pairing gap obtained from 
many-body calculations including medium polarisation effects and using realistic 
two- and three-body forces 0. 

For applications to nuclear masses, two phenomenological corrections are added 
to the HFB ground-state energy: (i) a Wigner energy and (ii) a rotational and vibra- 
tional spurious collective energy (see for instance Ref. |8| for details). The devi- 
ations between the 2149 measured masses of nuclei with N and Z > 8 given in 
the 2003 AME [4| and the predictions from our model HFB-17 are shown graph- 
ically in Fig. Q] The rms and mean values of these deviations are 0.581 MeV and 
-0.019 MeV, respectively. HFB-17 is the most accurate mass model ever achieved 
within the mean-field framework. The complete HFB-17 table of 8389 masses in- 
cluding all nuclei with Z,N > 8 and Z < 110 and lying between the proton and 
the neutron drip lines is available on our websiteQ. With the additional constraints 
on infinite nuclear matter, our HFB mass model is also particularly suitable for as- 
trophysics applications such as the description of neutron-star crusts Q. 

3 Applications to neutron stars 

Neutron stars are among the most compact objects in the Universe with a central 
density which can exceeds several times that found inside heavy atomic nuclei. Neu- 
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tron stars are born in the catastrophic gravitational core collapse of massive stars in 
supernova explosions. The outer layers of the star are formed of a solid crust ifTTIl . 
which at densities below the neutron drip threshold pnd — 4 x 10 11 g.cm~ 3 , is 
composed of a solid Coulomb lattice of neutron-rich nuclei with ZJA < 0.5 coex- 
isting with a degenerate gas of relativistic electrons. Our HFB mass models can be 
directly used to compute the composition and the equation of state of these layers 
following the classical work of Baym, Pethick and Sutherland [12]. We have found 
essentially the same sequence of nuclides (see for instance Table 4 of Ref. [2]) than 
that obtained by Haensel and Pichon lfT3l . In particular the last equilibrium nuclide 
at densities just below the neutron drip transition is 120 Sr. The inner crust of neutron 
stars, at densities above pnd up to about half the saturation density is permeated by 
a neutron ocean. The latter affects the properties of the "nuclei" by exerting a pres- 
sure on them and reducing their surface tension. Those nuclei are thus very different 
from those encountered on Earth. In order to make reliable predictions of the com- 
position of the inner crust, both nucleons bound inside "clusters" and free neutrons 
have to be described consistently. This can still be done using our effective force 
underlying our mass models by solving the HFB equations (0 with Bloch boundary 
conditions 

i/j[f (r +e,a)= cxp(ifc • e)ip[f (r, a) 

V^f (r + £, a) = exp(ifc • t)^f (r, a) (14) 

where k is the Bloch wave vector and £ is any lattice vector. But such calculations 
are computationally very expensive. So far self-consistent mean field calculations 
have been performed using a simpler approach based on the Wigner-Seitz (W-S) 
method with lfl4l and without pairing [15|. However the W-S treatment introduces 
spurious neutron shell effects which contaminate the results lfT6HT8l . For this rea- 
son, we have followed a different strategy by applying the Extended Thomas-Fermi 
method including proton shell corrections via the Strutinsky integral (see Ref. Q for 
details). This ETFSI method is not only a very fast approximation to Hartree-Fock 
equations, but it also avoids the pitfalls of boundary conditions than plagued current 
quantum calculations. This method could be similarly generalized to solve approx- 
imately the HFB equations (|5). But in our calculation of the equation of state H71 . 
pairing was neglected. Pairing is not expected to have any significant impact on the 
energy density and on the pressure because it only affects nucleon states lying close 
to the Fermi level. Indeed in uniform neutron matter, using our model HFB- 16 J8] 
we have found that the pairing contribution represents at most ~ 0.5% of the en- 
ergy per particle (without rest mass energy). But of course pairing is essential for 
studying neutron superfluidity. In the bottom layers of the crust where spatial inho- 
mogeneities are small, the effects of the nuclear clusters on the neutron superfluid 
can be estimated by solving the HFB equations perturbatively. Any field <fi(r) hav- 
ing the periodicity of the crystal lattice (i.e. single-particle fields, pairing field) can 
be expanded into Fourier series 



<t>(r) = (t>o + ^2(t>G exp(iG • r) 

G#0 



(15) 
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where G are reciprocal lattice vectors. The Fourier coefficients are defined by 

4>g = t^~ I d 3 r0(r)exp(-iG-r) (16) 

Vcell Jccll 

with Veen the volume of the unit cell. If <jf>(r) is spatially slowly varying, we will 
have \4>g\ *C |0o | f° r an Y G^O. Solving to lowest order the HFB equations (|9]l for 
neutrons with Bloch boundary conditions (fT4l i thus leads to the gap equation 
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where B n and v 71 " " are given by 

Bn = 7^-[ d 3 rB„(r), 3*n = -L- f dW[p n (r),p p (r)]. (18) 

Vcell Jccll ^ccll Jccll 

This equation is similar to the BCS Eq. < fT~3T > after substituting B n and n by their 
spatial average. This result is an illustration of the proximity effect: all particles 
whether inside clusters or not contribute to the pairing gap. Since the pairing gap is 
typically very small compared to the Fermi energy, we can approximate the chemi- 
cal potential by the latter 

A„ - B n kl n , k Fn = (37r 2 p„) 1/3 • (19) 

We have calculated the neutron pairing gap A n in the densest layers of the inner 
crust of neutron stars using the nucleon density profiles obtained with the ETFSI 
method. In order to study the modifications of the pairing gap due solely to the pres- 
ence of spatial inhomogeneities, we have applied our model HFB-16 [8 1 which was 
adjusted on the 1 So pairing gap of pure neutron matter, as calculated with realistic 
forces but without any medium effects. For comparison, we have also calculated the 
pairing gap A n of uniform neutron matter for the density p{ t corresponding to the 
density of free neutrons. Results are summarized in Table Q] We have found that 
nuclear clusters reduce the neutron pairing gap A n compared to A n . This can be 
understood by the fact that the neutron pairing gap arises from the spatial average 
of the pairing strength which is smaller inside clusters than outside. 



4 Conclusions 

Our HFB-17 mass model not only gives a better fit to the mass data than any 
other mean-field model, but is also by far the most microscopically founded. With 
the additional constraints on the properties of pure neutron matter, this model is thus 
expected to make more reliable predictions of highly neutron-rich nuclei. Besides it 
is very well suited for the description of astrophysical environments like supernova 
cores and neutron stars. In particular, this is the first of our models well adapted to 
the investigation of a possible superfluid phase in the inner crust of neutron stars. 
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Table 1. Equilibrium composition of the bottom layers of neutron-star crust and So neutron 
pairing gaps, p is the average nucleon density, Z and A the equilibrium numbers of protons 
and nucleons in the W-S cell respectively (as obtained from the ETFSI method [7| with 
effective force BSkl6 |8|), p„ the neutron density outside clusters, A n is the 1 Sq pairing gap 
of uniform neutron matter for the neutron density p£, A n is the neutron pairing gap obtained 
by solving the HFB equation perturbatively. 



P [fm" 3 ] 


Z 


A 


pi [fm~ 3 ] 


A n [MeV] 


A n [MeV] 


0.070 


40 


1258 


0.060 


1.79 


1.48 


0.065 


40 


1264 


0.056 


1.99 


1.72 


0.060 


40 


1260 


0.051 


2.20 


1.96 


0.055 


40 


1294 


0.047 


2.40 


2.21 


0.050 


40 


1304 


0.043 


2.59 


2.45 
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